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Fully-microscopic No-core Shell Model (NCSM) calculations of all stable s and p shell nuclei are used to 
determine a realistic NN interaction, JISP16, describing not only the two-nucleon data but the binding energies 
and spectra of nuclei with A < 16 as well. The JISP16 interaction, providing rapid convergence of the NCSM 
calculations, is obtained in an ab exitu approach by phase-equivalent transformations of the JISP6 NN interaction. 
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To complement the successful but computa¬ 
tionally intensive l ab initio ’ No-core Shell Model 
(NCSM) |T], we introduce the l ab exitu ’ NCSM. 
While the former has proven very successful 
for light nuclei when one includes three-body 
(NNN) forces |2I3| . the computational com¬ 
plexity motivates us to introduce an approach 
that simultaneously minimizes NNN forces while 
providing more rapid convergence with a pure 
nucleon-nucleon (NN) force. We invoke directly 
an end-goal of nuclear theory (hence the term ‘ ab 
exitu 1 ), a successful description of nuclear proper¬ 
ties, including the available NN data, to develop 
a new class of NN potentials that provide accu¬ 
rate descriptions of a broad range of nuclear data. 

To achieve this, we form a union of two re¬ 
cent techniques — the J-matrix inverse scatter¬ 
ing HSH and the NCSM [T], A major ingredient 
of our approach is the form of the NN interac¬ 
tion (a small matrix in the oscillator basis), which 
is chosen to provide rapid convergence of many- 
body observables within the NCSM. Indeed, we 
show below that results up through A = 16 ob¬ 
tained directly with the bare interaction (one that 
accurately describes the NN data) are close to 


those obtained with the effective interaction and 
are very useful to establish the confidence region 
for the binding energy. 

Since this is a departure from the more tra¬ 
ditional approach, we motivate our development 
with observations concerning the successful ab 
initio approaches to light nuclei. Indeed sev¬ 
eral promising microscopic approaches have been 
introduced and tested extensively with realistic 
NN interactions (see [7; and references therein) 
and with realistic NN + NNN interactions |8I21 
[31. Progress towards heavier nuclei appears lim¬ 
ited only by scientific manpower and by available 
computers. However, all approaches face the ex¬ 
ponentially rising computational complexity in¬ 
herent in the quantum many-body problem with 
increasing particle number and novel schemes are 
needed to minimize the computational burden 
without sacrificing realism and precision. 

The earliest and most successful in reaching 
nuclei beyond A = 4 is the Green’s-function 
Monte Carlo (GFMC) approach [8J whose power 
has been used to determine a sequence of ever- 
improving NNN interactions mum in conjunc¬ 
tion with highly precise NN interactions m that 
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fit a wide selection of low-lying properties of light 
nuclei up through A = 10. In addition, the 
Hamiltonians are tested for their predictions in 
infinite systems m- According to our usage of 
terminology here, the application of GFMC to 
determine successful NNN interactions is an ex¬ 
cellent example of an ab exitu approach. 

Now, we ask the question whether it is pos¬ 
sible to go even further and search through the 
residual freedoms of a realistic NN interaction to 
obtain new NN interactions that satisfy three cri¬ 
teria: (1) retain excellent descriptions of the NN 
data; (2) provide good fits to light nuclei; and (3) 
provide improved convergence properties within 
the NCSM. The challenge to satisfy this triad of 
conditions is daunting and we are able to provide 
only an initial demonstration at the present time. 

We are supported by the work of Polyzou and 
Glockle who demonstrated [T3] that a realistic 
NN interaction is equivalent at the A = 3 level to 
a realistic NN+NNN interaction where the new 
NN force is related to the initial one through a 
phase-equivalent transformation (PET). The net 
consequence is that properties of nuclei beyond 
A = 3 become dependent on the freedom within 
the transformations at the A = 3 level. It seems 
reasonable then to exploit this freedom and work 
to minimize the need for the explicit introduction 
of three and higher body forces. However, we do 
not surmise that we would be able to eliminate 
them completely. 

We start from the realistic charge-independent 
NN interaction JISP6 [6] that provides an ex¬ 
cellent description of the deuteron properties [6j 
and NN scattering data with y 2 /datum = 1.03 
for the 1992 np data base (2514 data), and 1.05 
for the 1999 np data base (3058 data) [14] . JISP6 
provides also a very good description of the spec¬ 
tra of p shell nuclei, but we find that it overbinds 
nuclei with A > 10. To eliminate this deficiency, 
we exploited PETs to modify the JISP6 in vari¬ 
ous partial waves. The resulting interaction, here¬ 
after referred to as JISP16 since it is fitted in our 
ab exitu approach to the spectra and bindings 
of stable A < 16 nuclei, can be obtained from 
the initial ISTP interaction in the same manner 
as JISP6 in Ref. [5] but with a different set of 
PET angles. These angles associated with uni¬ 


tary transformations (see Refs. [5|6j for details) 
mixing the lowest s and d oscillator basis states 
in the coupled sd waves and the lowest oscilla¬ 
tor basis states in the 3 p 2 , 3 pi, 3 Po, 3 ^2 and 1 pi 
waves are d = —11.0°, +5°, -6°, -10°, +25° and 
— 12° respectively. The JISP16 matrix elements 
in the oscillator basis with trio = 40 MeV that 
differ from those of JISP6, are presented in Ta- 
bles[l][2] The JISP16 predictions for the deuteron 
rms radius = 1.9643 fm and quadrupole mo¬ 
ment Q = 0.288585 fm 2 slightly differ from the 
JISP6 results since the JISP16 and JISP6 PET 
angle in the sd coupled waves is slightly differ¬ 
ent ($ = -11.0° versus —11.3°). In this paper, 
we include all NN partial waves up to l = 4 and 
include in Tables [T] [2] the added matrix elements. 

Our fitting procedure is one of ’trial-and-error’ 
where we worked with only a few partial waves 
that we deemed important for these nuclei. We 
fit only the excitation energies of few lowest 6 Li 
levels and the 6 Li and 16 O binding energies. To 
save time, we performed the NCSM calculations 
in small enough model spaces (up to 10 tko for 6 Li 
and up to 4 huj for 16 0). After obtaining a reason¬ 
able description of these observables, we checked 
that the binding energies and spectra of all the 
rest s and p shell nuclei are well-described in 
small model spaces. The results presented below 
are obtained in the ab initio NCSM calculations 
with the obtained NN interaction, the ab exitu 
JISP16, in larger model spaces. This description 
of the binding energies is somewhat worse than 
the one obtained during the fit in smaller model 
spaces, however it is still very reasonable. In a 
future effort, we will perform a thorough search 
through the space of possible PETs that should 
further improve the description of nuclear prop¬ 
erties while retaining the excellent description of 
the NN data. 

We illustrate our approach with the 16 O ground 
state energy in Fig. [TJ The variational princi¬ 
ple holds for the bare interaction results; hence 
the upper bound (UB) for the ground state en¬ 
ergy is the minimum of its hu> dependence in the 
bhu model space. In the calculations with the 
effective interaction obtained by the Lee-Suzuki 
transformation, the quoted result is convention¬ 
ally associated with the minimum of the Huj de- 
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Table 1 

JISP16 non-zero matrix elements in fun = 40 MeV units in the uncoupled NN partial waves that differ 
from the respective JISP 6 matrix elements and of the JISP16 matrices in higher partial waves. 


K 


i 

nn 


V 1 

v n,n-f- 1 y n+l,n 


v l = v l 

*n,n+2 v n+2,n 


V 1 
y > 


l 

nn 


V 1 _ L 1 =V’ i .i 

v n,n-f- 1 y n+l,n 


v l — v l 

v n,n +2 v re+2,n 


1 pi partial wave 


3 Po partial wave 


0 0.4864373541 -0.2359869829 0.3117643519 

1 -0.1487460250 -0.1438603014 

l g 4 partial wave 

0 -0.0159359974 0.0110169386 

1 -0.0229351778 0.0073206473 

2 -0.0056121168 


0 0.1571004930 -0.1425039101 0.2505691390 

1 -0.2172768679 -0.0981725471 

3 <?4 partial wave 

0 -0.0762338541 0.0498484441 

1 -0.1107702854 0.0371277135 

2 -0.0295683403 


Table 2 

Same as in Table Q] but for the coupled NN waves. 


sd coupled waves 



Vnn' matrix elements 


n 

y ss 

v nn 

Y ss —V ss 

v n,n+ 1 v n+l,n 


0 

-0.5125432769 

0.2139078754 



V dd , matrix elements 


n 

t rdd 
* nn 

v dd ,, = V d ii 

r n,n+1 r n+l,n 


0 

0.0551475852 

-0.0952367414 



V ad , = Ki'n matrix elements 


n 

v sd , =y ds , 

v n,n— 1 n—l,n 

t rsd _ T/'ds 

* nn 'nn 

v ad ,, =y d f n 

r n,n+l r n+l,n 

0 


-0.4035852241 

0.2003382771 

1 

-0.0464306332 



pf coupled waves 

VZ matrix elements 

n 

ypp 

v nn 

y p p =y«b 

r n,n+l r n+l,n 

t/PP _ yPP 

v n,n+ 2 * / n+2,n 

0 

-0.1933759934 

0.1508436490 

-0.1072949881 

1 

-0.0277262441 

0.0964883300 



Vnn' matrix elements 


n 

V pf ,=V fp , 

v n,n— 1 v n— l,n 

•S-e 

II 

y p/ ,, =y / P 1 

r n,n+1 r n+l,n 

0 


0.0195093232 

0.0020663826 

1 

-0.0252003957 

0.0236188613 



dg coupled waves 


V dd matrix elements 

n 

ydd 
* nn 

T/dd 

*n,n+l 

e 

? + 

II 

0 

-0.0226611102 

0.0231171026 

1 

-0.0514940563 

0.0256493733 

2 

-0.0329967376 

0.0061799968 

3 

-0.0002368252 




Vnn' matrix elements 


n 

V" 

v nn 

y 9 ® 

n,n+l 

— V" 

*n+l,n 

0 

0.0435654902 

-0.0276372780 

1 

0.0537629744 

-0.0140723375 

2 

0.0079901608 




V dp , matrix elements 


n 

ydg —ya d 

v n,n— 1 v n— l,n 

ydg . 

v nn 

is 

II 


0 

1 0.0791431969 

2 0.0660805779 

3 0.0029846726 


-0.0392683838 

-0.0874578184 

-0.0334474774 
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huj (MeV) 


Figure 1. (Color online) The fojj dependence of 
the 16 O ground state energy obtained with bare 
JISP16 and effective interaction based on JISP16 
in a sequence of N max huj model spaces up to 
Nmax = 6; the lines marked as Exp., UB and 
LB show the experimental ground state energy, 
the upper bound and the suggested lower bound 
for the NCSM ground state energy predictions. 


pendence. This minimum is seen from Fig. [T] to 
ascend with increasing model space. Based on our 
results in lighter systems with larger spaces that 
show uniform convergence of this minimum, the 
minimum obtained in the 6 Huj model space is a 
suggested lower bound (LB) for the ground state 
energy. The difference between these upper and 
lower bounds is our estimate for the ‘error bars’ 
of our predictions. These error bars suggest rea¬ 
sonable convergence is attained but this requires 
verification in larger basis spaces. 

Similar trends are found for most of the p shell 
nuclei. We present in Table [3] their binding ener¬ 
gies obtained with both bare and effective interac¬ 
tions. We also quote the hio values providing the 
minimum with the effective interaction. The dif¬ 
ference between the given result and the result ob¬ 
tained with the same Hlo in the next smaller model 
space is presented in parenthesis to give an esti¬ 
mate of the convergence of our calculations. We 
quote our differences in significant figures from 
the rightmost figure of the stated result, omit- 


Table 3 


Binding energies (in MeV) of nuclei obtained with 
bare JISP16 and effective interaction generated 
by JISP16. 


Nucleus 

Nature 

Bare 

Effective 

hco Model 

(MeV) space 

3 H 

8.482 

8.354 

8.496(20) 

7 

14 hu 

3 He 

7.718 

7.648 

7.797(17) 

7 

\vhu> 

4 He 

28.296 

28.297 28.374(57) 

10 

Whui 

6 He 

29.269 


28.32(28) 

17.5 

12fno 

6 Li 

31.995 


31.00(31) 

17.5 

Yihjj 

7 Li 

39.245 


37.59(30) 

17.5 

10 Hlo 

7 Be 

37.600 


35.91(29) 

17 

10 hco 

8 Be 

56.500 


53.40(10) 

15 

8 Huj 

9 Be 

58.165 

53.54 

54.63(26) 

16 

8Huj 

9 B 

56.314 

51.31 

52.53(20) 

16 

8Huj 

10 Be 

64.977 

60.55 

61.39(20) 

19 

8 Huj 

10 B 

64.751 

60.39 

60.95(20) 

20 

8 Huj 

10 C 

60.321 

55.26 

56.36(67) 

17 

8 Huj 

n B 

76.205 

69.2 

73.0(31) 

17 

6 Huj 

n C 

73.440 

66.1 

70.1(32) 

17 

6 Huj 

12 b 

79.575 

71.2 

75.9(48) 

15 

6hu> 

12 c 

92.162 

87.4 

91.0(49) 

17.5 

6Huj 

12 N 

74.041 

64.5 

70.2(48) 

15 

6 Huj 

13 b 

84.453 

73.5 

82.1(67) 

15 

6 Huj 

13 c 

97.108 

93.2 

96.4(59) 

19 

6 Huj 

13 N 

94.105 

89.7 

93.1(62) 

18 

6 Huj 

13 0 

75.558 

63.0 

72.9(62) 

14 

6 Huj 

14 c 

105.285 

101.5 

106.0(93) 

17.5 

6 Huj 

14 n 

104.659 

103.8 

106.8(77) 

20 

6 hu> 

14 0 

98.733 

93.7 

99.1(92) 

16 

6hu> 

15 N 

115.492 

114.4 

119.5(126) 

16 

6 Huj 

15 0 

111.956 

110.1 

115.8(126) 

16 

6 Huj 

16 0 

127.619 

126.2 

133.8(158) 

15 

6 Huj 
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Table 4 


Ground state energy E gs and excitation energies E x (in MeV), ground state point-proton rms radius r p 
(in fm) and quadrupole moment Q (in e ■ fm 2 ) of the 6 Li nucleus; Hu = 17.5 MeV. 


Interaction 

Method 

Nature 

JISP6 

NCSM, 10 Hu [6] 

JISP16 

NCSM, 12 Hu 

AV8’+TM’ 
NCSM, 6 Hu 0 

AV18+UIX 

AV18+IL2 

GFMC 8 15j 

GFMC HO 15 

£ ss (l^,0) 

-31.995 

-31.48 

-31.00 

-31.04 

-31.25(8) 

-32.0(1) 

r P 

2.32(3) 

2.083 

2.151 

2.054 

2.46(2) 

2.39(1) 

Q 

-0.082(2) 

-0.194 

-0.0646 

-0.025 

-0.33(18) 

-0.32(6) 

E x ( 3+0) 

2.186 

2.102 

2.529 

2.471 

2.8(1) 

2.2 

£*(0+1) 

3.563 

3.348 

3.701 

3.886 

3.94(23) 

3.4 

£*(2+0) 

4.312 

4.642 

5.001 

5.010 

4.0(1) 

4.2 

E x { 2+1) 

5.366 

5.820 

6.266 

6.482 


5.5 

£*(i£,o) 

5.65 

6.86 

6.573 

7.621 

5.1(1) 

5.6 


ting decimal points to save space. The ground 
state energy of A = 6, 7 and 8 nuclei converges 
uniformly from above with both the bare and ef¬ 
fective interactions. We present in Table [3] only 
the effective interaction results for these nuclei 
due to their superior convergence features. For 
these nuclei, an extrapolation based on the fit by 
a constant plus exponential function for different 
Hu values may be useful. For 6 Li, this extrap¬ 
olation results in a binding energy of 31.70(17) 
MeV where the value in parenthesis is the un¬ 
certainty of the fit. A similar extrapolation for 
6 He results in a binding energy of 28.89(17) MeV 
which is bound with respect to the a + n + n 
threshold. We note that the bare interaction re¬ 
sults for A = 6 nuclei are very close to the effec¬ 
tive interaction ones demonstrating a remarkable 
softness of the JISP16 interaction: the 6 Li and 
6 He binding energies are 30.94(44) and 28.23(41) 
MeV respectively, the extrapolations of the bare 
interaction bindings produce 31.33(12) MeV for 
6 Li and 28.61(12) MeV for 6 He. 

The nuclear Hamiltonian based on the ab exitu 
realistic NN interaction JISP16, is seen to re¬ 
produce well the binding energies of nuclei with 
A < 16. The lowest state of natural parity has the 
correct total angular momentum in each nucleus 
studied. The experimental binding energies of all 
nuclei presented in Table [3] either lie within error 
bars of our predictions or are close to our sug¬ 


gested LB based on the effective interaction cal¬ 
culations. Generally JISP16 slightly underbinds 
only nuclei in the middle of the p shell. The differ¬ 
ence between UB and LB is small, suggesting that 
JISP16 provides good convergence. However, our 
error bars increase as binding energy decreases in 
a chain of isobars (cf the results for 13 0 and 13 N). 

We present in Tables |4] and [5] spectra and 
ground state properties of 6 Li and 10 B which are 
known I2MT5ITTT1 to be sensitive to an explicit 
NNN interaction. Here, the ab exitu JISP16 
NN interaction alone provides a good descrip¬ 
tion. The JISP16 6 Li spectrum seems to be less 
favorable than that provided by our JISP6 in¬ 
teraction specifically fitted to the 6 Li spectrum. 
However, the JISP16 6 Li spectrum is competi¬ 
tive with those of realistic NN + NNN poten¬ 
tial models. Also, we obtain a good description 
of the 6 Li quadrupole moment Q that is a rec¬ 
ognized challenge due to a delicate cancellation 
between deuteron quadrupole moment and the d 
wave component of the a-d relative wave func¬ 
tion. We observe that Q and the point-proton 
rms radius r p have a more prominent Hu depen¬ 
dence than the binding energy. 

The 10 B properties are also seen to be well- 
described with the JISP16 interaction contrary 
to previous results from pure realistic NN inter¬ 
actions m- We note that the 10 B spectrum 
depends on Hu at N max = 8 but not so strongly 
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Table 5 

Same as in Table [4] but for the 10 B nucleus; Huj = 15 MeV. 


Interaction 

Method 

Nature 

JISP16 

NCSM, 8 Huj 

AV8’+TM’ 
NCSM, 4 Huj |2j 

AV18+IL2 
GFMC m 

Egs (3+, 0) 

-64.751 

-60.14 

-60.57 

-65.6(5) 

r P 

2.30(12) 

2.168 

2.168 

2.33(1) 

Q 

+8.472(56) 

6.484 

+5.682 

+9.5(2) 

£*(lf,0) 

0.718 

0.555 

0.340 

0.9 

E x ( 0+1) 

1.740 

1.202 

1.259 


£*(1^,0) 

2.154 

2.379 

1.216 


£*(2+0) 

3.587 

3.721 

2.775 

3.9 

£*(3+,0) 

4.774 

6.162 

5.971 


£*(2+1) 

5.164 

5.049 

5.182 


£*(2+0) 

5.92 

5.548 

3.987 


£*(4+,0) 

6.025 

5.775 

5.229 

5.6 

E x (2+,1) 

7.478 

7.776 

7.491 



as to alter our main conclusions. For example, the 
minimum of the 10 B ground state corresponds to 
Huj = 20 MeV while the minimum in the first ex¬ 
cited state energy occurs at Hu = 15 MeV. We 
present in Table [5] the 10 B properties obtained 
with Hco = 15 MeV, i.e. with the Hu> value cor¬ 
responding to the minimum of the first excited 
state since it has a more pronounced Huj depen¬ 
dence than the ground state. The 10 B ground 
state spin was not previously reproduced with a 
pure realistic NN interaction. We observe that 
our description of the 10 B spectrum is somewhat 
better than the one obtained with the Argonne 
AV8’ NN potential and Tucson-Melbourne TM’ 
NNN force. In particular, we reproduce the or¬ 
dering of 10 B levels except for the (3^,0) state. 
We note that the (3^,0) state is also too high 
with the AV8’ + TM’ interaction. 

In constructing ISTP J5], JISP6 [6] and JISP16 
potentials we adopted only the accepted symme¬ 
tries of the NN interaction and neglected explicit 
constraints such as the long-range behavior from 
meson-exchange theory. However, this does not 
mean that the JISP16 NN interaction is inconsis¬ 
tent with meson-theoretical forms of the NN in¬ 
teraction. On the contrary, it is well-known that 


the one-pion exchange (OPE) dominates the NN 
interaction in higher partial waves and the long- 
range behavior of NN interaction in lower par¬ 
tial waves. In this context, we showed in Ref. [5] 
that our scattering wave functions in higher par¬ 
tial waves are nearly indistinguishable from those 
of the Nijmegen-II OPE potential. Also, in lower 
partial waves, our wave functions are very close 
to those of Nijmegen-II at large distances and a 
small difference is seen only at higher energies. 
Finally, we introduced the PETs of JISP6 and 
JISP16 only in lower partial waves and only in 
a few lowest oscillator components of the poten¬ 
tial with a large value of Huj = 40 MeV. As a 
result, PETs reshape the wave functions at short 
distances (< 1 fm). Thus, the JISP16 interaction 
appears to be consistent with the well-established 
OPE tail as embodied in the Nijmegen-II NN in¬ 
teraction. 

We propose our ab exitu JISP16 as a realistic 
NN interaction since it describes the two-body 
observables with high precision. In addition, it 
provides a reasonable and economic description of 
properties of many-body nuclear systems in the 
microscopic NCSM approach. Economy arises 
from the softness of the interaction represented 
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in a separable oscillator form. Short distance 
phase-equivalent transformations adjust the off- 
shell properties successfully to reduce the roles of 
multi-nucleon interactions. The particular mech¬ 
anism of this reduction is not clear at the present 
time. However, our results as well as the success 
of the approach of Ref. EZ!> clearly demonstrate 
that such a mechanism exists and should be stud¬ 
ied in detail. We plan to study this with explicit 
NNN interactions. 

We conclude that the many-body nuclear 
Hamiltonian obtained in our ab exitu approach 
is realistic from the point of view of providing a 
good description of a wide range of nuclear data. 
The suggested JISP16 NN interaction opens a 
path for extending realistic microscopic theory to 
heavier nuclei, to achieve better convergence and 
to obtain improved agreement with experiment. 
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